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Abstract 



We evaluate the supersymmetry reach of the Large Hadron Collider within 
the gauge-mediated low energy supersymmetry breaking framework, assuming 
that a neutralino is the second lightest sparticle and that it decays promptly 
into a gravitino which escapes detection. We find that the maximum reach is 
obtained via a search for inclusive 77 + events coming dominantly from 
chargino and neutralino production. Assuming an integrated luminosity of 
10 we find that LHC experiments will be able to probe values of the 

model parameter A < 400 TeV, corresponding to nig < 2.8 TeV. A measure 
of the model parameter A may be possible from the photon spectrum. 
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Typeset using REVTgX 



Most phenomenological analyses |]] of supersymmetry (SUSY) have been carried out 
within the framework of the mSUGRA model ||[], where gravitational interactions induce 
effects of SUSY breaking (assumed to occur in a "hidden sector" ) into the "observable sector" 
of ordinary particles and their superpartners. In this case, the particle-sparticle mass gap 



in the observable sector is M Q b s ~ (F)/Mpi anc k, where J (F) is the scale at which SUSY 
is broken. To ensure that SUSY can stabilize the elementary scalar electroweak symmetry 
breaking sector, we require M ohs ~ 250 GeV, and so choose (F) ~ (10 11 GeV) 2 . The 
gravitino, which acquires a mass ~ M obs , then decouples from the observable sector, and is 
unimportant for particle physics phenomenology. 

Recently, however, there has been a resurgence of interest [P|-|i"2"| in the phenomenology 



of a class of models where usual gauge interactions, not gravity, serve as the messengers of 
SUSY breaking. Here, it is assumed that there is a new set of particles that directly couple 
to the SUSY breaking sector, and which also have non-vanishing SM interactions. The idea 
then is that these new particles, which first feel the effects of SUSY breaking, convey these 
down to the observable sector via their gauge interactions. This framework differs from the 
mSUGRA framework in two important ways. First, because SUSY breaking is fed down to 
the observable sector by gauge interactions, the mass difference between ordinary particles 
and their superpartners depends on their gauge interactions - gluinos are thus heavier than 
SU{2) gauginos which are heavier than binos, and rriq > m| ~ rriy > rni R . Second, since 
gravitational interactions are not the messengers of SUSY breaking, M obs is not suppressed 
by Mpi anc k as in mSUGRA. Instead, we now have sparticle masses ~ o(Fm) /M, where M is 
the mass scale of the messenger particles and (Fm) the induced SUSY breaking scale in the 
messenger sector, with a the relevant gauge coupling of the sparticle. The important thing 
is that it is possible to get weak scale sparticle masses with (Fm) ~ 100 TeV if M is 100- 
1000 TeV. Indeed it could well be that the intrinsic SUSY breaking scale (F) is of the same 
magnitude, and supersymmetry breaking is determined by 100 TeV scale physics. If this 
is the case, then the gravitino mass uiq = (F)/(v3Mp) (here M P = M P i anck /\^8n) could 
be 0(1 eV). This is important because (the longitudinal components of) such superlight 
gravitinos have phenomenologically relevant couplings to ordinary particles Jl3| . 



In the specific model that we adopt for our analysis, the messenger sector comprises of 
n 5 sets of "quark" and "lepton" superfields in a 5+5 representation of SU(5), with superpo- 
tential couplings to a gauge singlet superfield (S) which acquires vacuum expectation values 
for both its scalar (S) and auxiliary (Fm) components. The former results in masses for 
messenger quark and lepton superfields, whereas the latter is the result of the feed down of 
SUSY breaking from the hidden sector. Assuming that (Fm) -C (S) 2 , the gauginos acquire 
a mass 

m Xi = £^A, (1) 
whereas scalar components of chiral superfields get a mass, 



^scalar = 



(2) 



with A = (Fm) / (S), and a± given in terms of the usual hypercharge coupling g' by a± = 
Finally, C3 = | for colour triplets and zero for colour singlets while C2 = | for weak doublets 
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and zero for weak singlets. These relations, which are independent of the messenger sector 
superpotential couplings Aj get corrections of ~ (F M ) / (Xi(S 2 )) which are ignored in the 
subsequent analysis. The gravitino mass is given by, 

V3M P X(F M )V3M P - Gmv VSM p - [) 

The parameter Ccrav that we have introduced is > 1 since the SUSY breaking scale in the 
messenger sector cannot be larger than the intrinsic SUSY breaking scale in the hidden 
sector, and the superpotential coupling A < 1. SUSY breaking A-parameters and the B- 
parameter are induced only at two loops so that it is reasonable to suppose that these are 
small. The supersymmetric fi parameter is not determined by how SUSY breaking is medi- 
ated but is fixed (up to a sign) by radiative symmetry breaking as in the mSUGRA frame- 
work.Q A complete theory that includes the dynamics of SUSY breaking will presumably 
yield a value of /i consistent with this. Eq. (p and Eq. (fj) should be regarded as boundary 
conditions for the gaugino and scalar masses valid at the messenger scale M = X(S), so that 
these parameters need to be evolved down to the weak scale relevant for phenomenological 
analysis. Our framework is thus completely parametrized byf\ 

(A, M, n 5 , tan (3, sgn //, C Grav ) 

which together with Standard Model parameters completely determine the masses and cou- 
plings of all particles in the observable sector. The messenger mass scale M can range from 
A to Mqut or Mp. We focus our attention on models where M ~ 1 PeV, so that the SUSY 
breaking scale (Fm) (and possibly (F)) may be relatively low. 

In order to facilitate simulation of models with gauge mediated supersymmetry breaking 
(GMSB) we have recently included in ISAJET v7.37 a 'GMSB option' that allows one 
to use the GMSB parameter set introduced above as an input. ISAJET then evolves the 
gaugino and scalar masses as given by Eq. (|IJ) and Eq. (0) from the scale M down to the scale 
relevant for phenomenology, and computes the "MSSM parameters" that are then used in 
the evaluation of sparticle cross sections and decay patterns. The decays of neutralinos into 
gravitinos, Zi — > G7, Zi — > GZ and Zi — > Gh, H, A as well as (approximately) the Dalitz 
decay Zi — > e + e~G are included in ISAJET. The decays £r — > £G and fx — > tG, as well as the 
three body decays ]15j (mediated by a virtual neutralino) £p — > fifi and £p — > Tit£ have also 
been included. Widths for corresponding three body decays mediated by virtual chargino 
exchange are suppressed by the lepton Yukawa coupling, and have not yet been included. 



It is customary to eliminate tan/3 in favour of Bq, the B parameter at the scale M. Since Bq 
is small in this framework, this constrains tan (3 to be 20-30, the precise value depending on the 
choice of A. As in Ref. 0] , we do not impose this constraint since we expect that it will be altered 
by new interactions introduced to generate \x dynamically. 

2 It is possible to construct [Q] more complicated models which lead to different boundary conditions 
for scalar and gaugino masses from Eq. (|l]) and Eq. (^). We do not consider such scenarios in our 
analysis. 
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These can be significant only when — m^ x ^ m T so that the neutralino-mediated three 
body decays of /1r are kinematically very suppressed or forbidden. 

Within the GMSB framework, the parameter A sets the scale of sparticle masses. The 
dependence on the messenger scale (M) , which specifies where the boundary conditions ([!]) 
and (0) hold is presumably logarthmic. To help the reader link the model parameters with 
the more familiar sparticle masses, we show in Fig. 1 contours of rrig in the A — tan f3 plane 
for (a) n 5 = 1, and (6) n 5 = 2. We fix m t = 175 GeV, M = 1 PeV and take fj, > 0. We 
have checked that the spectrum is only weakly sensitive to M for a range of M between 
A and 1 PeV, and also that it is insensitive to the sign of /i except for very low values of 
tan/3. The SU{2) and U{1) gaugino masses are roughly given by ^rrig and \rrig, respectively. 
The scales of the sfermion masses can roughly be read off from Eq. @). In Region I below 
the diagonal solid line Q in frame ( a) , < , so that Z\ — > G7 (this dominates the 
decays to Z or Higgs bosons), so that isolated photons are a hallmark of SUSY within this 
framework. In Region II, sparticles cascade decay to fi (except immediately above this line 
where —m^ < m T ), which then decays to tG, so that for this portion of the parameter 
space SUSY events contain r leptons in addition to other leptons and jets 0. In frame 
(b), Zi is once again the next to lightest SUSY particle (NLSP) in Region I, while in the 
other regions f\ is the NLSP (though mg H and mp, R may be very close to m^). In Region 
II, m^ R > nig > m^, while in Regions III and IV, Z\ — > £r£ is kinematically possible. In 

Region III, £r — > £G because the decay £r — > rf\£ is kinematically inaccessible, while in 
Region IV, both these decays are possible and compete. The phenomenology will clearly be 
sensitive to just where we are in the parameter space. 

The grey region on top of each frame is excluded because proper electroweak symmetry 
breaking is not obtained: the wedge in the upper left corner is excluded because m| fl < 
while in the band on the top, m\ < 0. The various hatched regions are excluded by the 
non-observation of sparticles or a Higgs boson at LEP 2. In the region with vertical shading, 
< 60 GeV (in the corresponding portion of Region I, this limit may be stronger as the 
final state also contains photons); in that shaded horizontally, < 88 GeV; in the small 
diagonally shaded region at low tan/3, < 89 GeV (the current lower limit on the SM 
Higgs boson mass). While these limits \T7\ are essentially those obtained in the mSUGRA 
model, the strongest constraint comes from the non-observation of acollinear photon pairs 
from e + e~ — > Z X Z\ — > 77GG, and is special to this framework. The ALEPH collaboration 
P~8|] has excluded < 71 GeV assuming that Z\ ~ B and m iR = 1.5m^ , which is a good 



approximation for 715 = 1. This exclusion is shown as the diagonally shaded vertical strip in 



3 There are finite one loop corrections [16| which alter the relationship between the fermion mass 
and the corresponding Yukawa coupling that have not been included in ISAJET. These corrections, 
which can be important for large tan (5, would shift the exact location of this line. Thus, for model 
parameter values close to this line, these corrections could radically alter the phenomenology if the 
NLSP shifts from being a neutralino to being a stau. While these corrections are crucial to include 
in any study relating SUSY phenomena to the underlying parameters (that may be determined by 
physics at yet higher energy scales) of the theory, this issue is not critical to our main point which 
is a study of the LHC reach for the case where Z\ is the NLSP. 
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the figures.[] The DO experiment at the Tevatron from a non-observation of di-photon 
events with large B< T have inferred a lower bound ~ m yj> 1 > 150 GeV, which essentially 
excludes this same region. 

We have already noted that if model parameters happen to be in Region I of Fig. 1, in 
addition to multijet and multilepton events expected in the mSUGRA framework, SUSY 
events would contain, in addition, hard isolated photons, provided that the SUSY breaking 
scale (F) is not much larger than (Fm) for the messenger sector, i.e., that C grav is not too 
large. In Region II (and III and IV in frame b), SUSY events would contain additional r 
(and possibly e or fi) instead of photons. The main purpose of this study is to assess the 
SUSY reach of the Large Hadron Collider within the GMSB framework, assuming that the 
parameters are in Region I. For this reason, we will mainly focus on the case = 1 in the 
rest of the analysis and defer the study of the phenomenology where sleptons are lighter 
than Z\ to a subsequent analysis. 

We begin by showing the production cross sections as a function of the parameter A for 
a variety of sparticle production reactions in Fig. 2. We have fixed M = 1 PeV, n 5 = 1, 
> which we take to be our default values, and taken tan/5 = 2. The cross sections 
are independent of Ccrav We use CTEQ3L structure functions |2(J for our analysis. We 
see that squark and gluino production dominates for small values of A < 175 TeV, beyond 
which electroweak production of (lighter) charginos and neutralinos is the dominant source 
of sparticles. The production of charginos or neutralinos in association with gluinos or 
squarks is always sub-dominant. The dotted curve, which includes cross sections for ££, £v 
and vv production, shows that slepton pair production is an order of magnitude smaller 
than chargino-neutralino production. Except when A ~ 50 — 100 TeV where v£l production 
is also important, slepton production is dominated by £r£r production {£ = e,/z in the 
figure), a reflection of the fact that 2m ^ ~ ~ mp. It is interesting to observe that £r£r 
production which results in the essentially SM background-free ££77 + I$ T events may be 
observable (if they can be separated from other SUSY sources of dileptons) for A values as 
large as 225 TeV, corresponding to = 400 GeV. 

For detector simulation at the LUC, we use the toy calorimeter simulation package IS- 
APLT. We simulate calorimetry covering —5 < 77 < 5 with cell size Arj x A<p = 0.05 x 0.05. 
We take the hadronic energy resolution to be 50%/\^E © 0.03 for \rj\ < 3, where © denotes 
addition in quadrature, and to be 100%/^/E © 0.07 for 3 < \rf\ < 5, to model the effective 
Pt resolution of the forward calorimeter including the effects of shower spreading, which is 
otherwise neglected. We take electromagnetic resolution to be 10%/y/E © 0.01. Although 
we have included these resolutions, which are typical of ATLAS and CMS f2"2fl , we have 
made no attempt to estimate the effects of cracks, edges, and other problem regions. Much 
more detailed detector simulations are needed to understand the effects of such regions and 
of the resulting non-Gaussian tails, particularly on the ]$ T resolution. 

Jets are found using fixed cones of size R = y/ Ar/ 2 + A<p' 2 = 0.7 using the ISAJET routine 
GET JET. Clusters with Et > 100 GeV and 1 77 (jet ) | < 3 are labeled as jets. Muons and 



4 For 715 = 2 in Fig. lb, although we have shown the same limit of 71 GeV, the reader should be 
aware that this is overly conservative because, for any given value of m~ , mg R will be lighter than 
for n.5 = 1, and the Z1Z1 cross section will be correspondingly larger. 
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electrons are classified as isolated if they have pr > 10 GeV, \r)(£)\ < 2.5, and the visible 
activity within a cone of R = 0.3 about the lepton direction is less than Et (cone) = 5 GeV. 
Photons with > 10 GeV and \r]\ < 2.5 are considered to be isolated if there is less than 
2.5 GeV of other activity in a cone of R = 0.3 about the photon. 

To enhance the SUSY signal over SM background processes, we require that each event 
contains at least two entities (an entity is a jet, an electron, a muon or a photon) with 
E T > 100 GeV and I$ T > 100 GeV and at least one isolated photon with E T > 40 GeV. The 
inclusive two photon sample is required to have a second isolated photon with E T > 25 GeV. 
As in previous analyses we find it useful to use variable cuts for the two hardest 

entities and I$ T : we require E?i, E?2, If? > E^, where we vary the parameter E^. 

To estimate SM background we have not done a new simulation, but used our previous 
computations of backgrounds for SUSY signals without photons to obtain an estimate of the 
background level for the present situation. The background cross sections as obtained from 
Fig. 3 of Ref. are shown in the first three rows of Table 1 for several values of E^. In 
this computation, the basic lepton and jet identification requirements were identical to the 
ones used in the present analysis and Ex(jl), Et(]2) and I$ T were all required to be larger 
than E^. Radiation of each hard isolated (i.e. large angle) photon is suppressed by a factor 
a (we neglect coefficients of 0(1) and factors of 2ir), so that we may estimate the physics 
backgrounds for SUSY events with one or two isolated photons by multiplying the total 
background in Ref. []24j] by a and a 2 , respectively. This background level is shown in the 
next two rows of the Table. What about non-physics backgrounds from mis-identification 
of jets or leptons as photons? We estimate this by assuming that the jet gamma rejection 
is 1/5000 |25| , and that SUSY events typically have a jet multiplicity ~ 10. The chance of 
an electron to fake a photon is much smaller (~ 10 -5 ) so that this background is negligible. 
We may thus obtain a crude estimate of the non-physics backgrounds to single (two) photon 



events by reducing the total background in Ref. |24j] by a factor of 500 (tM^ ~ 2 x 10 6 ) 



This is shown in the last two rows of Table 1. An important conclusion that one draws from 
the Table is that while the detector-dependent background to single photon events may be 
comparable to the physics background, it is relatively unimportant for the inclusive two 
photon SUSY event sample. 

The reader may worry that the background sample always contains two jets with Et > 
Ej,, while a portion of the signal may instead contain leptons or photons instead of jets as 



the primary entities, the background to which has not been included in Ref. yA . While 



we have not estimated this, we do not expect it to be a problem since it is reasonable to 
suppose that in the SM, events with very hard isolated leptons and/or photons are much 
rarer than those with hard jets. 

In order to estimate the SUSY reach of the LHC, we have computed the signal cross 
sections for the one and two-photon inclusive SUSY samples. Typically, we found that 
after acceptance cuts, the two photon inclusive signal rate is suppressed relative to the 
corresponding single photon rate by a factor of 2 — 3, depending on sparticle masses. Since 
two photon backgrounds are much more strongly suppressed the best reach is obtained in 
the 77 + Iftrp + X channel. We show LHC cross sections for such events versus A after 
cuts in Fig. 3 for (a) jj, > 0,tan/9 = 2, (b) jj, > 0,tan/3 = 35, and (c) /i < 0,tan/5 = 2. 
The parameters and M are set at our default values. In our computation, we assume 
that the decay Z\ — > 7G occurs instantaneously; i.e. the photon emerges from the primary 
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vertex. (The typical decay length is millimeters for C grav = 1. This is short enough to 
be unimportant for photons but is easily measured via Z\ — > e + e~G Dalitz decays.) Also 
shown is the background level, estimated as described above. We show cross sections for 
E T = 100 GeV (solid) and 200 GeV (dashed). The scale on the top axis is the gluino mass, 
which is essentially the same for each of the frames. We see that by choosing the larger value 
of E T , especially for A > 200 — 250 GeV, the SM background can be greatly reduced with 
only modest loss of signal. Indeed, the signal is then only rate limited, and we obtain a 5 
event level, which we take to define the reach of the LHC, with an integrated luminosity of 
10 fb~ l if A < 400 TeV (corresponding to rrig < 2.8 TeV) in all three frames. The expected 
SM background is 0.5 event, so that the Poisson probability of fluctuation to give more than 
five events is 1.4 x 10~ 5 . Of course, one must view this number in perspective, considering 
the crudeness of our background estimate. Note however, that even if the background is 
larger than our estimate by an order of magnitude, the typical experiment would exclude 
the SM at 99% CL if A < 380 TeV, or a "5a discovery" would be possible for A < 300 TeV. 

We remark that we have also checked for several cases along the boundary of Region I 
in Fig. lb that the SUSY signal is observable even if n.5 = 2, as long as Z\ — > G7 well inside 
the detector. We expect that the signal is observable over this entire region. 

Finally, we briefly consider prospects for determining the parameter A that sets the 
sparticle mass scale in our scenario. At first, it might seem that the mean values of the total 
scalar energy, the visible invariant mass, or the E? of the hardest entities in SUSY events 
would scale with A. This is, however, not the case because the dominant SUSY sources 
change from g and q to (much lighter) charginos and neutralinos as A increases from below 
200 TeV to above 200 TeV. The mass of Z\ (which is dominantly a bino), on the other hand 
always scales with A, so that it is reasonable to expect that the px of the photons produced 
via its decays to the gravitinos would scale similarly. To illustrate this we show in Fig. 4 
the pt distribution of the hard photon in SUSY events where both photons are detected for 
several values of A. For definiteness, we have illustrated this for four values of A with other 
parameters as in Fig. 3 a, and for Ej, = 100 GeV. Assuming that our background estimate is 
valid, we see that the distributions show a clear scaling with A, though for the highest value 
of A an integrated luminosity of 30-50 fb^ 1 would be required to construct this distribution. 
The pt distribution of the second photon shows the same scaling behaviour. 

Although the Pt{i) distributions in Fig. 4 scale nicely with A, it is not enough to de- 
termine A as 715 is not known. The most obvious strategy for the determination of is 
to obtain the ratio of sfermion and gaugino masses. This may, for example, be possible by 
focussing on event chains in which the decays Z 2 — > ££r — > IIZ\ — > ££jG (for n 5 = 1) 
or even Z\ — > £r£ — > UG (for n 5 = 2) can be isolated. Additional information may also 
be available from the overall mass and Et of the hardest entities in SUSY events. The 
distribution of decay lengths of the NLSP, which may be possible to determine from the 
displacement of the secondary vertex can yield information about the underlying SUSY 
breaking scale. A study of these issues, and also the reach of the LHC when the NLSP is 
not the lightest neutralino, is in progress. 
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TABLES 



TABLE I. SM background cross sections from Fig. 3 in Ref. \24\ as a function of Eff, defined 
in the text. The 21 background includes contributions from both same sign and opposite sign 
leptons. Trilepton background cross sections are much smaller. The next two rows, which are 
obtained by multiplying the total background by a ~ 10~ 2 and a 2 ~ 10 -4 , give our estimate of 
physics sources for single photon and two photon inclusive events passing the SUSY selection cuts. 
Detector-dependent backgrounds, which are estimated as discussed in the text, are shown in the 
last two rows. 



E C T (GeV) 


100 


200 


300 


a(E T ) (fb) 


10 4 


500 


50 


a(U) (fb) 


1200 


30 


4 


a(2£) (fb) 


700 


10 


0.4 


a x a Tot (fb) 


100 


5 


0.5 


a? x a To t (fb) 


1 


0.05 


0.005 


Fake(l 7 ) (fb) 


20 


1 


0.1 


Fake(2 7 ) (fb) 


0.02 


10" 3 


10" 4 
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FIG. 1. Various regions in the A — tan/3 plane for M = 1000 TeV and fi > for a) = 1 
and 6) ri5 = 2. The grey region is excluded by theoretical considerations while the various hatched 
regions are excluded by experimental constraints discussed in the text. Also shown by the dotted 
curves are contours of gluino mass. In Region I, Z\ is the NLSP while in the other regions, the 
NLSP is a slepton. 



12 



10 
10 
10 
10 
10 

1 

10 



-1 





n 5 =1 




gg + gq + qq 




fi > 








tan/? =2 




II 








Assoc. Prod. 


1 1 1 1 1 


'*''*** ^""N. 

I I I I I I I I I I 


I'r-k 


sJ '"'r-t-j.. i i i i i i 



50 100 150 200 250 300 350 400 450 500 

A (TeV) 

FIG. 2. The production cross sections versus A for various sets of SUSY processes at a 14 TeV 
pp collider within the GMLESB framework. We show the cross sections for M = 1000 TeV, with 
other parameters as shown in the figure. For sleptons, we sum over the first two families. 
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FIG. 3. Signal cross sections for inclusive 77 + If T + X events at the LHC after the cuts 
described in the text for two values of the cut parameter Ej.. We fix M = 1000 TeV. Our estimate 
of the background level from SM sources is shown by the lines denoted by bg. The five event level 
is also shown for an integrated luminosity of 10 fb^ 1 . 
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FIG. 4. The transverse momentum distributions for the hard photon in inclusive 77 + I$ T 
events at the LHC as computed within the GMSB framework for several values of A and the cut 
parameter E%. = 100 GeV. We choose M = 1000 TeV. 
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